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We present large-scale density-functional theory (DFT) calculations and temperature programmed desorption 
measurements to characterize the structural, energetic and vibrational properties of the functionalized molec- 
ular switch 3,3',5,5'-tetra-fert-butyl-azobenzene (TBA) adsorbed at Au(lll). Particular emphasis is placed 
on exploring the accuracy of the semi-empirical dispersion correction approach to semi-local DFT (DFT-D) in 
accounting for the substantial van der Waals component in the surface chemical bond. In line with previous find- 
ings for benzene and pure azobenzene at coinage metal surfaces, DFT-D significantly overbinds the molecule, 
but seems to yield an accurate adsorption geometry as far as can be judged from the experimental data. Com- 
paring the trans adsorption geometry of TBA and azobenzene at Au(l 1 1) reveals a remarkable insensitivity of 
the structural and vibrational properties of the — N=N— moiety. This questions the established view of the role 
of the bulky tert-butyl-spacer groups for the switching of TBA in terms of a mere geometric decoupling of the 
photochemically active diazo-bridge from the gold substrate. 



I. INTRODUCTION 

The ultimate goal of nanotechnology is control over 
molecular-scale mechanical and electronic components. An 
oft-proposed route to this goal is the construction of com- 
ponents from controllable single molecules. An important 
class of such molecules with obvious applications is formed 
by those that have properties that are reversibly and bi-stably 
modifiable by external stimuli - so-called molecular switches. 
One example is the azobenzene molecule (HjCg-lSHSI-CeHs), 
which can be bi-stably photo-isomerized between its planar, 
C2/1 symmetric trans and torsioned-twisted, C2 cis conform- 
ers in both solution and gas-phase. The high yield and sta- 
bility of this reaction have rendered azobenzene an archetype 
of molecular switch research with proposed technical appli- 
cations inclu ding e.g. light-driven actuator^ and information 
storage media^. 

For many such applications, switching of molecules at solid 
interfaces - adsorbed at metal surfaces, for example - is 
of particular interest. Unfortunately however, the switching 
properties of azobenzene have proven highly sensitive to the 
adsorbate-substrate interaction: Even at nearly chemically in- 
ert close-packed noble metal surfaces, switching of surface- 
adsorbed azobenzene by light has never been achieved, and 
switching by excitation with a scanning tunneling microscope 
(STM) tip has been successful only at Au(lll'/. A natural 
route to restoring the adsorbate switching ability is to fur- 
ther decouple the frontier n,n and K* orbitals responsible^-^ 
for the gas- or liquid phase photo-isomerization from the 
substrate electronic structure. With these frontier orbitals 
largely located at the central diazo (— N=N— ) bridge an intu- 
itive idea to achieve such a decoupling is to functionalize the 
molecule with bulky spacer groups that prevent a closer en- 
counter of the photochemically active unit with the substrate. 
This is precisely the notion behind the arguably to date most 
studied such adsorbate, 3,3',5,5'-tetra-ferf-butyl-azobenzene 
(TBA)^"-^. TBA consists of azobenzene functionalized with 
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FIG. I: Perspective views of trans and cis azobenzene and TBA, 
together with an illustration of the diazo-bridge bond length i^nn and 
the dihedral CNNC angle CO. The latter is defined as the smallest an- 
gle between two planes spanned by the — N=N— bridge and — C— N— 
and — N— C— bonds to the first and second phenyl-ring, respectively. 
The C atoms of the phenyl-rings have been darkened to allow easier 
distinction of azobenzene backbone and functional butyl groups in 
TBA. 



four ferf-butyl (-C-(CH3)3) groups in the phenyl-ring meta po- 
sitions as illustrated in Fig. [T] These 'table legs' were indeed 
found to enhance the switching efficiency of the adsorbed 
species, as e.g. indicated by the successful TBA switching 
by light at Au(lll)''3nil. However, attempts at STM-tip in- 
duced switching at ostensibly comparable substrates such as 
Ag(lll)'^ and Au(lOO)^^ have been unsuccessful, indicating 
that TBA is not significantly more robust to specifics of the 
substrate interaction than pure azobenzene. 

These circumstances beg the question, how and to what de- 
gree the TBA butyl groups really 'decouple' the photochromic 
moiety from the substrate. For a corresponding atomic-scale 
understanding the detailed characterization of the adsorbate 
geometry and binding constitutes a prerequisite and is the 
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main objective of the present contribution. In contempo- 
rary surface science, corresponding analyses are increasingly 
performed by quantitative first-principles electronic structure 
calculations. Particularly density-functional theory (DFT) 
with present-day local or semi-local exchange-correlation (xc) 
functionals has developed into an unparalleled workhorse for 
this task, with often surprising accuracy particularly with re- 
spect to structural properties of the surface adsorption system. 
For TBA at Au(lll) corresponding calculations are already 
challenged by the extension of the functionalized molecule 
and the simultaneous necessity to describe the metal band 
structure within a periodic supercell approach. On a more fun- 
damental level, the real limitation comes nevertheless from 
sizable dispersive van der Waals (vdW) contributions to the 
surface chemical bond as characteristic for organic molecules 
containing highly polarizable aromatic ring systems^'*. With 
local and semi-local xc functionals inherently unable to ac- 
count for such contributions'^ skewed, if not qualitatively 
wrong results must therefore be suspected. We have recently 
quantified thi s for azobenzene at the close-packed coinage 
metal 

surface JSEl. 

Comparison to detailed structural and 
energetic reference data from normal-incidence x-ray stand- 
ing wave (NI-XSW) and temperature programmed desorp- 
tion (TPD) measurements for azobenzene at Ag(l 1 1) demon- 
strates that the prevalent semi-local DFT xc treatment leads 
indeed to a significant underbinding with key structural pa- 
rameters deviating by more than 0.5 A'^. 

For system sizes as those implied by the adsorption 
of azobenzene or TBA an appealing and computationally 
tractable possibility to improve on this situation is a semi- 
empirical account of dispersive interactions within the frame- 
work of so-called DFT-D schemes . In this approach the 
vdW interactions not described by present-day xc function- 
als are approximately considered by adding a pairwise inter- 
atomic CeR^^ term to the DFT energy. At distances below 
a cut-off, motivated by the vdW radii of the atom pair, this 
long-range dispersion contribution is heuristically reduced to 
zero by multiplication with a short-range damping function. 
While the applicability of this approach to adsorption at metal 
surfaces is uncertain {vide infra), our recent benchmark study 
for azobenzene at Ag(lll) revealed that in particular the 
most recent DFT-D scheme due to Tkatchenko and Scheffler 
(TS)^^ yields excellent structural properties, albeit at a notable 
overbindin^iS. 

In this contribution we further explore the generality of 
this finding by analyzing the adsorption geometry, vibra- 
tions and energetics of TBA at Au(lll). Comparison 
against our recent near edge x-ray absorption fine structure 
(NEXAFS)^^ and high-resolution electron energy loss spec- 
troscopy (HREELS)^^ measurements, as well as a complete 
analysis of new TPD data confirms the accurate structural and 
vibrational predictions reached by the DFT-D TS scheme. The 
again obtained significant overbinding furthermore supports 
the interpretation'** that the neglect of metallic screening is 
the main limitation in the application of this scheme to the ad- 
sorption of organic molecules at metal surfaces. Comparing 
the TBA data to those for pure azobenzene at Au(l 1 1) we find 
a qualitatively different adsorption geometry for the functional 



backbone in the case of the cis isomer. For the trans isomer, 
on the other hand, we determine an intriguing structural and 
vibrational insensitivity of the photochemically active central 
diazo-bridge to the presence of the bulky spacer groups. The 
role of the latter for the switching efficiency is therefore more 
subtle than simple geometric decoupling and will be the topic 
of continuing work in our groups. 



II. METHOD 
A. Theory 

The DFT-D methodology followed in this work corre- 
sponds exactly to that employed in our preceding work for the 
adsorption of azobenzene. We therefore restrict ourselves here 
to a brief account and refer to our previous publications'^ '^ 
for an in-depth description of the underlying concepts and 
technical details. 

The DFT calculations were performed using a plane- wave 
basis set and library ultrasoft pseudopotentialP^ as imple- 
mented in and distributed with the CASTEP^^ code. The gen- 
eralized gradient approximation (GGA) to the xc-functional 
with the parametrization suggested by Perdew, Burke and 
Ernzerhof^^ (PBE) was used throughout. In the spmt of the 
DFT-D approach, the lack of vdW interactions in this semi- 
local functional is approximately corrected with an additional 
analytical, two-body inter-atomic potential. At long range, 
this potential equals the leading Cg^^^-term of the London 
series, where R is the inter-atomic distance, and Cg a so-called 
dispersion coefficient. At short range, this long range poten- 
tial is matched to the DFT inter-atomic potential by a damp- 
ing function f{R,R'^), typically modulated by the vdW radii 
R^ of the atom pair In this work we use the material-specific 
C(, and /?" parameters, as well as the damping function form 
suggested by Tkatchenko and Scheffler^^ (henceforth denom- 
inating corresponding results as PBEh-TS). This scheme ac- 
counts to some degree for the bonding environment through 
a Hirshfeld-analysis based adjustment of the Ce parameters, 
which in our previous work on azobenzene at coinage metal 
surface ga ve a s uperior performance compared to other DFT- 
D scheme^EMl. 

The analytical form of the dispersion correction poten- 
tial brings the advantage that dispersion-corrected total ener- 
gies of geometries that are fully relaxed with respect to the 
dispersion-corrected forces can be obtained (employing an in- 
house extension to the CASTEP code) at essentially the same 
computational cost as a regular GGA calculation. On the other 
hand the assumption of a simple two-body form for the dis- 
persion potential inherently neglects the effect of electronic 
screening of the vdW interaction^^, which particularly for the 
here studied adsorption at metal surfaces is expected to lead 
to an overestimation of the binding energy '^. A second limi- 
tation of the semi-empirical DFT-D approach might arise for 
adsorbate molecules which also interact covalently with the 
substrate. This typically leads to molecule-substrate bond dis- 
tances that are so short that the uncertainties in the heuristic 
damping function of the dispersion term may mingle in an un- 
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controlled way with deficiencies of the employed semi-local 
DFT xc functional. 

The surface calculations were performed within supercells, 
using (11 l)-oriented metal slabs with (6x5) surface unit-cells 
and at least 18 A vacuum. We verified that lateral interactions 
between adsorbed TBA and its periodic images are negligble 
at the GGA-PBE level. In the dispersion correction poten- 
tial corresponding interactions were also switched off, so that 
our calculations should give a fair account of TBA adsorption 
in the low-coverage limit. As in our preceding work we ne- 
glected the subtle effects of the long-range Au(lll) surface 
reconstruction. Full geometry optimizations (to within resid- 
ual forces below 30 meV/A) of all molecular degrees of free- 
dom were correspondingly performed with TBA adsorbed on 
one side of a four layer bulk-truncated slab. Test calculations 
with appropriately saturated ferf-butyl groups indicated only 
a weakly expressed site specificity of these TBA functional 
groups at Au(lll). This suggests that the photochemically 
active diazo-bridge plays a prominent role in anchoring the 
molecule at the metal substrate. The geometry optimizations 
were correspondingly initiated with the TBA diazo-bridge lat- 
erally placed as in the previously determined optimal adsorp- 
tion geometry of azobenzene, corresponding to a 1:1 N-metal 
atom coordination^^. Harmonic vibrational spectra of the ad- 
sorbed species in these relaxed geometries (and in the gas- 
phase) were subsequently obtained from numerical Hessians 
calculated by finite differences and neglecting any degrees of 
freedom of substrate atoms. To efficiently parallelize over the 
432 displacements required for each spectrum, we have inter- 
faced with the Atomistic Simulation Environment including 
the 'phonopy' extension to iPSI. Not aiming to reproduce the 
HREELS intensities'^, simple Lorentzian broadening with a 
width of 1 meV was applied to the spectra for better visual- 
ization. 

For the energetic and electronic structure calculations the 
thus determined relaxed adsorbate geometries were inverted 
in the bottom layer, forming inversion-symmetric seven layer 
slabs with adsorbates at both sides. This zeroes the internal 
dipole moment of the slab and results in a substantially im- 
proved substrate electronic structure. The two central ener- 
getic quantities obtained with the resulting structures are the 
adsorption energy 

£'ads = 2 [^azo@(lll) ^-£'(111)] ^-£'azo(gas) ) (1) 

and the relative stability of adsorbed cis (C) and trans (T) con- 
former 

A^C-T = 2 ['^azo@(lll)(C) -£'azo@(lll)(T)] . (2) 

Here £'azo@(iii) is the total energy of the relaxed, double- 
sided azobenzene-surface system, £(iii) the total energy of 
the clean slab, and £'azo{gas) the total energy of the corre- 
spondingly relaxed gas-phase isomer (computed within F- 
point sampled (35 A x 45 A x 35 A) supercells). Where ap- 
plicable, TS DFT-D corrections calculated in optimized four 
layer slab and gas-phase geometries are added to £azo@(iii) 



and £'azo(gas)' respectively. Both quantities were also con- 
sistently zero-point energy corrected with the previously ob- 
tained vibrational frequencies. In the convention of Eq. ([T]) 
the adsorption energy of either cis or trans isomer at the sur- 
face is thus measured relative to its stability in the gas-phase 
at both pure PBE and dispersion corrected PBEh-TS levels of 
theory, and a negative sign indicates that adsorption is exother- 
mic. Consistently, a negative sign of AEq-t indicates a higher 
stability of the cis isomer Convergence tests show that at 
the employed plane wave cutoff of 450 eV and (2 x 3 x 1) 
Monkhorst-Pack (MP) grid^^ both energetic quantities are nu- 
merically converged to within ±30 meV. 



B. Experiment 

The TPD measurements were carried out under ultrahigh 
vacuum conditions. The Au(lll) crystal was mounted on 
a liquid nitrogen cooled cryostat, which in conjunction with 
resistive heating enables temperature control from 90 K to 
750 K. The crystal was cleaned by a standard procedure of 
Ar+ sputtering and annealing. The TBA was dosed by means 
of a home-built effusion cell held at 380 K at a crystal temper- 
ature of 260 K. In the TPD experiments, the substrate was re- 
sistively heated with a linear heating rate of 1 K/s and desorb- 
ing TBA was detected with a quadrupole mass spectrometer at 
the TBA-fragment mass of 190 amu (3,5-di-ferf-butyl-phenyl 
ion). This procedure was repeated for different dosing times 
corresponding to different initial TBA-coverages. 

As further discussed below, three desorption features (tti- 
as) are observed in the TPD. They are assigned to desorp- 
tion from the multilayer (ai) and the first monolayer (ML) 
(a2 + 0C3), where a2 represents the desorption of ss 10% of 
the monolayer coverage (for details see Fig. |5]and Ref. [TOll. 
The NEXAFS and HREELS measurements of Refs. I24I25I 
discussed in Section III were performed at a coverage of 1.0 
and 0.9 ML, respectively, which were prepared by heating the 
multilayer-covered surface to 340 K or to 420 K. 



III. RESULTS AND DISCUSSION 

A. Adsorption Geometry 

Our previous work on azobenzene at coinage metal surfaces 
points to an understanding of the surface chemical bond in 
terms of a balance of four major contributions: A covalent 
bond between the diazo-bridge and the metal, the vdW attrac- 
tion between the metal and the phenyl-rings, the Pauli repul- 
sion between the phenyl-rings and the metal, and the ener- 
getic pe nalty d ue to the distortion of the gas-phase molecular 
geometryl^El Yhis understanding should largely carry over 
to TBA at Au(lll), with the bulky tert-butyl groups partic- 
ularly adding to the vdW attraction and the molecular defor- 
mation upon adsorption. Key structural parameters to char- 
acterize the adsorption geometry for both trans and cis iso- 
mer are therefore the location and orientation of the central 
diazo-bridge, as well as the orientation of the planes of the two 
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FIG. 2: Side view of adsorbed cis TBA at Au(l 11), illustrating key 
structural parameters defining the adsorption geometry (see text): 
The vertical height Zi of the two diazo-bridge N atoms, as well as 
the out-of -horizontal phenyl plane bend angles 
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TABLE I: Structural parameters as defined in Figs. [T]and|2] Val- 
ues for adsorbed TBA at PBE+TS level of theory are compared to 
corresponding data obtained at PBE level of theory, as well as for 
adsorbed pure azobenzenei^. Additionally shown are the out-of- 
horizontal phenyl plane bend angles CO/ as determined recently by 
NEXAFS measurementJSI. 



phenyl-rings with respect to the surface. As further illustrated 
in Figs. [T]and|2]we concentrate on the — N=N— bond length 
t/NN and the vertical N atom— surface plane distances Zi for 
the prior. Here, the indices / — 1 and / = 2 represent the value 
for the N atom connected to the phenyl-ring closer to and fur- 
ther away from the surface, respectively. Obviously, zi ~ zi 
reflects an diazo-bridge that is oriented parallel to the surface, 
as we obtain for the more symmetric trans isomer through- 
out. To specify the position of the phenyl-rings we focus on 
the CNNC dihedral angle co and the out-of-horizontal bend 
angles co, of the two ohenyl planes, with the same convention 
for the index i— 1 , 2 to distinguish the closer and more distant 
phenyl-ring in the asymmetric adsorption geometry of. Fig. 
|2] Adapting to the convention employed in Refs. |f8|and|24l a 
dihedral angle of CO = 180° indicates a planar TBA molecule, 
while <bj ~ 0° denotes a phenyl-ring that lies parallel to the 
surface plane. 

Table |l] summarizes the detailed geometric parameters ob- 
tained from our calculations with the exception of the diazo- 
bridge bond length. For the latter we consistently compute 
only insignificant changes away from the TBA gas-phase 
value of 1.29 A (trans) and 1.28 A (cis) at both PBE and 
PBE+TS level of theory. While this indicates an overall mi- 
nor activation of the — N=N— bond upon adsorption, it also 
demonstrates an insensivity to the degree of vdW interactions 
accounted for in the calculations. Starting the more detailed 
presentation with the trans isomer, Fig.[3]displays a perspec- 
tive view of the overall TBA adsorption geometry. In line 
with the interpretation of STM^ ii surface vibrationaP^and 
NEXAFS^"* data, the central azobenzene moiety essentially 
maintains its gas-phase planarity, with its long axis aligned to 
about 5° along the direction of close-packed atom rows on the 




FIG. 3: Perspective view of the trans TBA adsorption geometry at 
Au(lll). 



(Ill) surface. Also consistent with the preferential orienta- 
tion deduced from NEXAFS the tert -butyl groups contact the 
surface with one C-H bond pointing towards the substrate. At 
the pure PBE level of theory, the central diazo-bridge coiTe- 
spondingly floats at a rather high height of about 4 A parallel 
to the surface. This is further away than for pure azobenzene, 
where this height was 3.5 A^. With the PBE functional not 
providing attractive vdW components to the molecule-surface 
interaction, the main effect of the bulky tert-butyl groups con- 
forms therefore with the intuitive expectation to simply lift the 
functional azobenzene backbone further away from the sur- 
face. 

This picture is obviously prone to change, once some ac- 
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FIG. 4: Perspective view of the cis TBA adsorption geometry at 
Au(lll). 



count of vdW attraction is added to the theoretical description. 
At the PBE+TS level of theory, the height of the diazo-bridge 
is indeed significantly reduced to 3.22 A, which is neverthe- 
less still large on the scale of a typical N-transition metal bond 
length. In contrast, the differential interaction with the surface 
is not significantly altered by the additional bonding compo- 
nent. The slight bending of the azobenzene moiety away from 
the ideal gas-phase planarity is basically the same in both 
the PBE and PBEh-TS adsorption geometry. Here, the par- 
allel orientation of the diazo-bridge to the surface plane and 
the corresponding slight upward bending of both phenyl -rings 
by ~ 10° are in very good agreement to the recent NEXAFS 
experimentJSI Indirectly, this corroborates some of the as- 
sumptions made in the NEXAFS data analysis, without which 
twice as large tilt angles would have resulted^'*. On the other 
hand, due the insensitivity of both structural parameters due 
to the degree of dispersive interaction in the calculations this 
agreement with experiment does unfortunately not allow any 
conclusion as to the accuracy of the description provided by 
the semi-empirical TS scheme. This is particularly unfortu- 
nate, as the latter predicts an intriguing similarity of the diazo- 
bridge height for TBA (3.22 A) and pure azobenzene (3.28 A), 
cf. Table [l] If correct, this obviously largely contradicts the 
afore mentioned intuitive view of the bulky spacer groups in 
terms of 'table legs' that decouple the photochemically active 
unit in a geometric sense. 

The similarity between functionalized and pure molecule 
does not extend to the adsorption geometry of the cis isomer 
at Au(lll). As shown in Fig. |4]this geometry is skewed for 
TBA with both the position and orientation of the two phenyl- 
rings distinctly different. This feature is consistently obtained 
at both PBE and PBEh-TS level of theory, while for pure 
azobenzene both theoretical descriptions agreed on an adsorp- 
tion structure with only a small sideward tilt away from a C2 
rotational symmetry around the central diazo-bridge. In the 
asymmetric adsorption mode of TBA the lower phenyl-ring 
is tilted out of the surface plane, yet without significant tor- 



sion that would place its two butyl groups at different heights 
above the surface. This is much different for the upper phenyl- 
ring, which does not only stand essentially upright, but is so 
torsioned that one of its butyl groups points largely towards 
the diazo-bridge and the other one away from it. Overall, the 
internal structure of the azobenzene backbone in this adsorbed 
state is therefore very similar to that in gas-phase TBA with 
even the CNNC dihedral angle co not much affected by the sur- 
face potential. Another noteworthy feature of the skewed ad- 
sorption mode is the prominent tilt of the central diazo-bridge, 
i.e. in contrast to adsorbed trans and cis azobenzene and ad- 
sorbed trans TBA the two N atoms are at a notably different 
vertical height, cf. Table [l] Qualitatively, this overall struc- 
ture is perfectl y consi stent with the understanding reached in 
the recent STlvF™^ and NEXAFS^"* measurements. As to 
the NEXAFS, the agreement is even quantitative in all re- 
spects. The determined out-of-surface-plane bend angles of 
both phenyl-rings agree very well with the assignments made, 
cf. Table U This also extends to the inclination of the cen- 
tral CNNC plane with respect to the surface plane^^. For this, 
NEXAFS determines about 60°, close to the 61° and 68° de- 
termined at the PBE and PBEh-TS level, respectively. 

To some extent unfortunate and similar to the situation for 
the trans conformer, none of these qualitative features, as well 
as tilt and inclination angles are very sensitive to the descrip- 
tion of vdW interactions in the calculations. Both PBE and 
PBE+TS give essentially identical results. The major fea- 
ture introduced again by the account of vdW attraction is an 
essentially rigid downward shift of the entire molecule by 
about 0.9 A. This is slightly larger than is the case for trans 
TBA 0.8 A) and leads at least for the lower N atom of the 
diazo-bridge to a vertical height above the sur face (2.34 A) 
that is close to the value determined for pure c/.s-azobenzene 
(2.23 A). With such prominent differences between PBE and 
PBE+TS, measurements of the vertical heights as was done 
by NI-XSW for pure azobenzene'^ would therefore again pro- 
vide a critical benchmark to judge on the accuracy of the ac- 
count of vdW interactions introduced by the semi-empirical 
TS scheme. In turn, however, precisely the lacking sensitiv- 
ity of the qualitative features, as well as tilt and inclination 
angles of both trans and cis adsorption geometries suggests 
that they are not much affected by the shortcoming of the em- 
ployed semi-local xc functional with respect to dispersive in- 
teractions. This increases the confidence that a quite reliable 
understanding of the adsorption geometry of both conformers 
has been reached by the present calculations. 



B. Energetics 

In view of the preceding study on azobenzen J^l it is clear 
that the binding energetics provided by the semi-empirical 
DFT-D approach deserves particular scrutiny. We performed 
TPD measurements to obtain an experimental reference value 
for the binding energy of the thermodynamic ally favored trans 
TBA on Au(lll). Figure |5] shows the TPD spectra as a func- 
tion of TBA surface coverage recorded at the fragment mass of 
190 amu ((C4H9)2-C6H^), with a linear heating rate of 1 K/s. 
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FIG. 5: Thermal desorption spectra (raw data) of TBA adsorbed on 
Au(ll 1) at different coverages as recorded with a linear heating rate 
of 1 K/s at the fragment-mass of 190 amu ((C4H9)2-C6H4"). The 
inset shows the desorption rate (ln{dries/dt) plotted against the re- 
ciprocal temperature at a coverage 6 of 0.1 monolayer (ML). The 
slope of the line, -E^^^/R with R the gas constant, then determines 
the activation energy for desorption isjes ^t this coverage, in this case 
£des(0.1 ML) = 1.65±0.08 eV. 



In the low coverage regime a broad desorption peak (as) is ob- 
served around 542 K which extends to lower temperature with 
increasing coverage, as has also been observed for TBA on 
Ag(ll 1 F^an d other azobenzene (derivatives) on noble metal 
surfacesSm xhis behavior is attributed to repulsive inter- 
actions between the adsorbed molecules (for example due to 
dipole-dipole interactions). 

After saturation of the desorption peak as a second feature 
(X2 develops at 394 K. Further increase in coverage leads to 
saturation of this peak and to the appearance of a new feature 
around 308 K labeled as ai. The ai peak increases in height 
and width with increasing coverage, showing a typical zero- 
order desorption behavior (data not shown). We therefore as- 
sign the ai peak to desorption from the multilayer, while 0.2 
and as are associated with desorption from the monolayer 
The a.2 component represents the desorption of about 10 % of 
the monolayer coverage and is attributed to desorption out of 
a densely packed TBA structurJ'^. 

In order to derive the activation energy for desorption, fides- 
of TBA on Au(l 11) in the low -coverage regime from this TPD 
data we utilize the so-called complete analysis^^, which has 
the advantage that knowledge about the absolute coverage is 
not required. For this analysis a family of TPD curves are 
measured as a function of initial coverage (0,) as shown in Fig. 
[5] These curves are subsequently used to construct a family 
of 9(r) curves via °^ rii^^^{t)dt = 9 with the peak area 
and rdes the desorption rate. Due to the known linear heat- 
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TABLE II: Energetics of TBA in gas-phase, and adsorbed at 
Au(l 1 1). All numbers are in eV. 



ing rate (in our case 1 K/s), this corresponds to a knowledge 
of 9(7^5) as a function of 9,-, where is the surface tempera- 
ture. Following this, an arbitrary coverage value 9i is chosen 
that is contained in each of the desorption curves. The des- 
orption rate at this coverage, rdeslQi), and the temperature at 
which this rate was obtained, Ti , are then read off from each 
of these desorption curves. Plotting ln[riies(9i)] versus l/Ti 
(Arrhenius-plot) finally allows to determine fides (9 1) as fol- 
lows directly from the Polanyi-Wigner equation^^. This is ex- 
emplarily shown for a coverage of 0.1 ML in the inset of Fig. 
|5] which yields an activation energy of 1.65 ±0.08 eV. 

In the low-coverage regime < 0.2 ML this analysis de- 
termines an essentially constant desorption energy of fides — 
1.70±0.1eV. Assuming no additional activation barrier for 
adsorption, this then provides the aspired benchmark num- 
ber against which the calculated fiads for trans TBA may be 
referenced. As apparent from Table [H] the values obtained at 
the PBE and PBEh-TS level are very consistent with the find- 
ings of the previous stud ies of benzene and pure azobenzene 
at coinage metal surface^^^^^ The lack of vdW interactions 
in the semi-local PBE functional leads to an overall negligi- 
ble binding of both TBA isomers, with the actual values for 
fiads in fact essentially identical to those computed for pure 
azobenzene at Au(l 11) before'^. Adding dispersive attraction 
within the DFT-D approach these binding energies are dramat- 
ically increased, such that in the end the intended dispersion 
'correction' amounts to more than 90% of the total binding 
energy. Compared to the experimental reference the PBEh-TS 
approach overbinds the trans isomer by about as much as 
PBE underbinds. In absolute numbers this is a disconcert- 
ing deviation of more than 1 eV. In our previous work we had 
assigned the corresponding overbinding determined for pure 
azobenzene to the neglect of metallic screening in the DFT-D 
approach'**. This argument was largely based on the intrigu- 
ing accuracy of the PBEh-TS adsorption geometry compared 
to the bond distances derived from NI-XSW measurements. 
As discussed above, such a conclusion on the determined TBA 
adsorption geometry is presently not possible, as all hitherto 
measured structural parameters are not very sensitive to the 
additional vdW attraction provided by the PBEh-TS scheme. 
On the other hand, the additional butyl groups lead to a signif- 
icantly increased overall binding energy of TBA compared to 
azobenzene, and in both cases PBE+TS overbinds by roughly 
40%(azo: 1.71 eV vs. 1.08 eV TBA: 3.00 eV vs. 1.70 eV). 
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Such a rather geometry-unspecific deviation between theoret- 
ical and experimental data quite well fits the anticipated effect 
of an overestimated uniform background potential. The latter 
is attributed to the neglect of screening of vdW contributions 
from more distant substrate atoms. 

For gas-phase TBA PBEh-TS yields a somewhat smaller 
cis-trans energy difference A/sc-t than for pure azobenzene, 
0.29 eV compared to 0.49 eV respectively. This arises pre- 
dominantly from the additional vdW attractions due to the 
butyl groups in the bend cis conformer, cf. Fig. [T] and is 
not found at the PBE level of theory, where A£c-t — 0.58 eV 
for both TBA and azobenzene. Compared to the gas-phase 
reference, the stability difference of the two TBA conformers 
is with 0.99 eV substantially increased at Au(l 1 1) at PBEh-TS 
level of theory. This is primarily due to the larger vdW at- 
traction possible in the planar trans adsorption mode and was 
equivalently obtained for pure azobenzene adsorptiorPl. This 
interpretation in terms of vdW is corroborated by the essen- 
tially unchanged AZic-x of gas-phase and adsorbed TBA at 
PBE level of theory, where these bonding contributions are 
absent, cf. Table |ll] Not withstanding, in light of the sus- 
pected overestimation of the vdW attraction within PBEh-TS 
we also expect that this increase of A/ic-T upon adsorption 
is overestimated. Extrapolating the roughly 40 % overshoot 
seen in the trans TBA binding energy, we would therefore 
cautiously conclude on only a moderate ^ 0.3 — 0.4eV in- 
crease of the cis-trans energy difference upon adsorption of 
TBA at Au(lll). While the therewith implied higher stabil- 
ity of adsorbed trans TBA is consistent with the existing ex- 
perimental data, there are to date unfortunately no dedicated 
measurements of A£c-t against which this estimate could be 
compared. 



C. Vibrations 

The picture arising from the geometric and energetic char- 
acterization of TBA at Au(lll) points to a surface chemical 
bond that is predominantly due to unspecific dispersive inter- 
actions. Particularly for the trans conformer, the structural 
similarity of the azobenzene and TBA adsorption complex 
with respect to the photochemically active diazo-bridge moi- 
ety (together with the in both cases negligible — N=N— bond 
activation) questions the oft-quoted function of the tert-butyl 
'spacers' in terms of a geometric decoupling from the surface. 
Complementary information confirming this view can come 
from an analysis of the surface vibrational modes. Figure [6] 
compares corresponding results computed for trans TBA at 
Au(lll) at PBE+TS level of theory with experimental data 
from HREELS^^. At first glance, we find overall agreement 
with all major features nicely reproduced by the calculations. 
As had already been noticed in the experimental work, the vi- 
brational spectrum exhibits only minor changes between ad- 
sorbe d cis a nd trans TBA and even TBA in the condensed 
pjj^gg, 25135136] rpjjjg similarly obtained in the calculations, as 
illustrated by the spectrum for cis TBA at Au(l 11) also shown 
in Fig. |6] 

The detailed inspection of the C-H vibrational modes in 



Trans. PBE 




Cis, Exp. 



50 100 150 200 250 300 350 400 
Energy (meV) 

FIG. 6: Comparison of calculated surface vibrational modes of trans 
(top, red lines) and cis (center, blue lines) TBA at Au(lll) at PBE 
(dashed lines) and PBE+TS (solid lines) level of theory to the corre- 
sponding HREELS spectra from ref. |25] (bottom, black solid lines. 
Upper: trans. Lower: cis). Center inset: The highest energy peaks 
of the spectrum, shown at larger scale. Note the C-H stretch peak on 
the left. The spectra have been vertically displaced for clarity. Not 
aiming to reproduce the HREELS intensities, the computed modes 
are convoluted with a Lorentzian broadening of 1 (meV). 



the energy range of 350 to 400 meV reveals nevertheless the 
subtle influence of the Au(lll) substrate on the vibrational 
spectra of the adsorbed species. Based on the calculations the 
three frequency bands centered at 368, 378, and 390 meV are 
assigned to the symmetric and asymmetric C-H vibrations of 
the tert-butyl groups, as well as the C-H bending modes at 
the phenyl-rings, respectively. There is a clear change of the 
former two vibrational bands when going from the PBE to 
the PBE+TS description, cf. the inset of Fig. |6] The overall 
reduced distance of the TBA molecule to the gold substrate 
upon inclusion of the attractive vdW forces leads to a clear 
splitting of the C-H vibrations of the fert-butyl moieties. In 
the adsorbed species the frequencies of the C-H bonds point- 
ing towards the substrate are red shifted by 5-lOmeV and 
the corresponding vibrations develop low-energy sidebands. 
These subtle but important changes are confirmed by the ex- 
perimental HREELS spectra. Due to the mainly in-plane char- 
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acter of the C-H phenyl-ring bending modes we observe only 
a small shoulder at 390 me V. In contrast the symmetric and 
asymmetric C-H stretch vibrations of the ferf-butyl legs at 368 
and 380 meV result in pronounced dipole-active bands. Upon 
trans to cis isomerization the number of C-H bonds point- 
ing towards and noticeably contacting the substrate is roughly 
halved, cf. Fig. [2] This leads in both experiment and theory 
to a reduction of the low-energy side bands best seen for the 
symmetric C-H stretch vibrations which constitute the lowest 
band at 368 meV. As the splitting of the C-H stretch vibrations 
happens only at a reduced distance to the substrate (PBE vs. 
PBEh-TS adsorption structure) it confirms the importance of 
vdW interaction in the bonding of TBA. Unfortunately this 
is not a quantitative benchmark for the accuracy of our semi- 
empirical PBEh-TS approach. 

More insight into the specific bonding of the diazo-bridge 
moiety can specifically come from a detailed analysis of the 
— N=N— stretch mode. However, identically obtained with 
PBE and PBEh-TS, the computed value of 184meV in ad- 
sorbed trans and cis TBA again reflects a negligible activation 
of the NN bond at the surface. Compared to the respective 
gas-phase values, the mode red-shifts only by about 3 meV 
for both isomers. Correspondingly, the obtained shift between 
this stretch mode in trans and cis TBA at Au(l 1 1) is very sim- 
ilar to the one computed for free TBA, 6 vs. 5 meV, respec- 
tively. This in turn compares very well to the 8meV mea- 
sured for T BA in the condensed phase with IR and Raman 
spectroscopjEHm. Such a negligible change of the cis-trans 
stretch frequency difference upon adsorption is not found for 
pure azobenzene at Au(lll). Here, the closer encounter of 
the diazo-bridge in the cis adsorption geometry softens the 
stretch mode by more than 6 meV. With the mode unaffected 
in adsorbed trans azobenzene, this leads to a concomitant in- 
crease of the cis-trans stretch frequency shift compared to the 
gas-phase. In cis TBA at Au(lll) a corresponding softening 
does not occur as the diazo-bridge does not come as close to 
the surface in the skewed adsorption mode. In this respect, 
one can really attest some geometric decoupling effect to the 
butyl groups, yet only for the cis isomer. For the trans iso- 
mer though, an equivalent analysis of the diazo-bridge against 
surface stretch or other torsional azo modes shows always the 
same similarity between TBA and azobenzene at Au(lll) as 
for the just discussed — N=N— stretch mode. Also the charac- 
terization of the surface vibrational properties supports there- 
fore the understanding that the effect of the bulky ferf-butyl 
groups for the switching properties of trans TBA at Au(lll) 
must be more subtle than a mere geometric decoupling of the 
central photochemically active molecular moiety. 



IV. SUMMARY AND CONCLUSIONS 

A methodological motivation of this work was to further 
explore the capabilities of the semi-empirical dispersion cor- 
rection approach to semi-local DFT in describing the adsorp- 
tion of complex organic molecules at metal surfaces. For this 
we have presented a detailed characterization of the geomet- 
ric, energetic and vibrational properties of the trans and cis 



isomers of the azobenzene derivate TBA at Au(lll). The 
findings are in all respects in line with the experience from 
preceding studies on benzene and pure azobenzene at coinage 
metal surfaces'^" '**: The additional account of attractive vdW 
interaction introduced by the PBE+TS scheme leads to a sig- 
nificant modification of the adsorption geometry, primarily in 
terms of bringing the molecule closer to the surface. The con- 
comitantly dramatically increased binding energy is notably 
overestimated compared to the reference value from TPD 
measurements. For azobenzene at Ag(lll) a corresponding 
overbinding - in fact in relative terms very much comparable 
to the one found here for TBA at Au(lll) - was attributed 
to the neglect of metallic screening of dispersive interactions 
in the semi-empirical DFT-D approachP^. This argument was 
largely based on the very accurate PBE+TS adsorption ge- 
ometry as compared to detailed structural data from NI-XSW 
measurements. Such a conclusion is, unfortunately, not yet 
possible in full for the here studied TBA at Au(lll). The 
qualitative features of the determined trans and cis adsorption 
geometries and even detailed tilt and incli nation angles agree 
all very well with existing data from STM ^I" ' ^^ ! HREELSSI 
and NEXAFS^^ experiments. The energetic lowering of the 
C-H stretch vibrations of the ferf-butyl groups pointing to- 
wards the substrate is only observed in the PBE+TS approach. 
This points towards the importance of the vdW component 
introduced by the DFT-D approach. However, no experi- 
mental reference data exist to date for the vertical height of 
the adsorbed molecule, that would allow to directly confirm 
the present working hypothesis, namely that the PBE+TS ap- 
proach is a very useful and computationally tractable tool to 
provide accurate structural information for adsorbed complex 
organic molecules. 



If this working hypothesis proves correct and the deter- 
mined PBE+TS adsorption geometries are indeed accurate, 
the prevailing preconception of the role of the bulky func- 
tional groups for the experimentally observed improved iso- 
merization ability of adsorbed TBA needs revision. Moti- 
vated by the gas-phase structure, this view discusses the butyl 
groups as spacers that help to decouple the molecular switch 
from the metal substrate. For the TBA cis isomer we indeed 
compute a skewed adsorption mode, in which one phenyl- 
moiety is much more tilted, up to the point of standing es- 
sentially perpendicular to the surface plane. This is quite dif- 
ferent to the more symmetric adsorption mode of cis azoben- 
zene at Au(lll), such that here one can attest some geomet- 
ric decoupling from the surface due to the butyl 'table legs'. 
Quite distinctly, the almost planar adsorption mode of trans 
TBA is very much comparable to the adsorption geometry of 
pure azobenzene. Particularly for flie photochemically active 
diazo-bridge moiety this similarity goes even down to essen- 
tially unchanged structural and vibrational properties. In this 
respect, a core message arising from the present study is that 
especially for the photo-switching of trans TBA the effect of 
the bulky spacer groups must be more subtle than the antici- 
pated mere geometric decoupling of the functional backbone 
from the metal surface. 
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